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ABSTRACT 

We present and analyze the orbits of eight nearby planetary nebulae (PNs) using two different 
Galactic models. The errors of the derived orbital parameters are determined with a Monte Carlo 
method. Based on the derived orbital parameters, we find that Sh 2-216, DeHt 5, NGC 7293, A21, 
and Ton 320 belong to the thin-disk population, and PC 1034-1-001 and A31 belong to the thick-disk 
population. PuWe 1 probably belongs to the thick-disk population, but its population classification 
is very uncertain due to the large errors of its derived orbital parameters. The PN-ISM interactions 
are observed for the eight PNs in our sample. The position angles of the proper motions of the PNs 
are consistent with the directions of the PN-ISM interaction regions. The kinematic ages of PNs are 
much smaller than the time for them to cross the Galactic plane. Using the models of Borkowski et 
al. and Soker et al., the PN-ISM interaction can be used to derive the local density of ISM in the 
vicinity of evolved PNs. According to the three-dimensional hydrodynamic simulations of Wareing et 
al. (WZO), Sh 2-216, A21, and Ton 320 are in the WZO 3 stage, PC 1034-hOOl and NGC 7293 are 
in the WZO 1 stage, and PuWe 1 is in the WZO 2 stage. 

Subject headings: planetary nebulae: general — planetary nebulae: individual (Sh 2-216, PuWe 1, 
A21, Ton 320, A31, PC 1034-^001, NGC 7293, DeHt 5) — star: white dwarfs 



1. INTRODUCTION 

Planetary nebulae (PNs) are the descendants of 
asymptotic giant branch (AGB) stars. According to the 
Interacting Stellar Wind (ISW) model (|Kwok et al.lll978l : 
lKwokl[l"982l ) , an AGB star loses most of its envelope mass 
with a slow stellar wind (AGB wind) until the hot core 
is exposed. After reaching the PN nucleus stage, a new 
faster wind (post- AGB wind) sweeps up the circumstel- 
lar material originating in the slow wind and forms the 
high-density and ionized nebular shell. Most PNs are 
not spherical and exhibit complex shapes: round, ellip- 
tical, butterfly, bilobed, and irregular. Stellar rotation, 
magnetic field , and binarity are important ingredien ts in 
shaping PNs ()De Marcoll2"00llBalick fc Frank|[2002l ). 

The interaction of a moving PN with t he interstel- 
lar rn edi um (ISM) was fi rst d iscussed by iGurzadvanl 
()1969( ). iBorkowsk i et all (| 19901 . hereafter BSS) exam- 
ined the morphology of nine PNs with a center star 
(CS) whose proper motions exceeded 0.015 arcsec yr~^ 
and found that seven of the examined PNs interact with 
the ISM. They also examined the morphology of large, 
nearby PNs, and f ound that many were asymmetric. 
iTweedv fc Kwitted (fT996h presented a CCD imaging at- 
las of large PNs (> 8') and found that 21 of 27 PNs 
in their sample revealed significant interaction with the 
ISM. Using deep narr owband CCD i mages of the nine 
largest, nearby PNs, iXilouris et al.l (|1996[ ) found that 
eight PNs in their sample revealed features indicating an 
interaction between the PNs and the ISM. More recently, 
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iRamos-Larios fc Phillips! (|2009l ) presented mid-infrared 
images for five PNs which appear to interact with the 
ISM. Usin g the dat a of the I NT Photometric Ha Survey 
(IPHAS) (|Drew et a l. 2005), ISabin et al.l ([20loh identi- 
fied 21 objects as possible PNs interacting with the ISM. 

On the theoretical side, the first theoretical st udy of 
the in teraction between PNs and ISM was made bv lSmitbl 
|1976|). Using a thin-shell ap proximation and simple an- 
alytic formulae, iSmithI ([1976') found that the motion of a 
PN through the ISM would be expected to cause a shift of 
the CS toward the leading surface of the interacting PN 
and ISM. BSS used simple analytical estimates to study 
the interaction between a moving PN and the ISM and 
found that the interaction produces asymmetries in the 
surface brightness of the PN. Using a two - dimen sional 
numerical hydrodynamic code, iSoker et al.l (|1991l here- 
after SBS) simulated the interaction between the PN 
and the ISM. They considered two limiting cases: adi- 
abatic and isothermal models. In their adiabatic model, 
a PN moves with large relative velocity (> 100 km s~^) 
within the ambient ISM with low density so that the 
cooling time is larger than the fiow time. The adia- 
batic model is relevant for PNs moving in the Galac- 
tic halo. In the isothermal model, which is relevant for 
PNs in the Galactic disk, the PN moves through the 
high-density ISM with moderate relative velocity ~ 60 
km s^^ and the cooling is very effective. The isothermal 
model of SBS valida tes the t hin-sh ell approximation orig- 
inally developed bv iSmithI ()1976D . SBS also found that 
the simple analytic model breaks down in their adiabatic 
model in which a Rayleigh- Taylor (RT) instability devel- 
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ops, leading to shell fragmentation, while the isothermal 
model shows no sign of the RT instability. iSoker fc Dganil 
(|1997f) studied the interaction of a PN with a magnetized 
ISM and showed that the ISM magnetic field, which was 
not incorporated in the numerical simulations of SBS, 
can lead to the development of RT instability even in 
the isothermal environment. 

According to the simple thin-shell model for the 
PNs moving through the ISM, the interaction be- 
tween PNs and ISM can be divided into three phases 
(jBorkowski et a l. 1990; Sokcr ct al. 1991). First, when 
the density of the young PN is much larger than that 
of the ambient ISM, the PN will expand freely. Sec- 
ond, when the PN density drops below a critical value, 
the ISM pressure becomes comparable to the one in the 
PN shell and the upstream side of the PN shell will be 
compressed and the surface brightness in the interaction 
region will increase. When the density in the PN shell 
has dropped to a lower value, the shape of the PN shell 
will be significantly distorted and the CS of the PN will 
move close to the upstream edge. Finally, the CS of the 
PN will move outside of the PN shell, and the PN shell 
will be stripped and mixed with the ISM. 

In the PN-ISM interaction models of BSS and SBS, 
they have assumed that the PN shell has formed be- 
fore the interaction with the ISM and the interaction 
only be observable during the late stage of the PN evo- 
lution. Moreover, high ISM density, large relative ve- 
locity, or both explain the observed a symrn etry in the 
external shell of a PN. iVillaver et al.l (|2003f l performed 
two-dimensional numerical simulations of the evolution 
of an AGB star moving through its surrounding ISM. 
They found that during the early stage of the evolution 
of an AGB star, the slow AGB wind can interact with 
the ISM. The PN forms later when the fast post-AGB 
wind sweeps up the stellar material in the slow AGB 
wind. A spherical innermost shell shaped by the fast 
post-AGB wind and an external s hell with a bow shock 
character can form. In simulations IVillaver et al.l (|2003[ ) 
adopted two low ISM densities and a very conservative 
relative velocity of the PN, 20 km s"^ Thus, high ISM 
density, large relative velocity, and a magnetic field are 
not necessary to explain the observed asymmetries in the 
external shells of PNs. Moreover, interaction can occur 
even during early AGB evolution, which is in agreement 
with the observed common asymmetries relating to the 
interaction with ISM in most of the halos of PNs, not 
just in those of old PNs. 

Considering a wide range of P N velocities, mass-loss 
rates of the CS and ISM densities, iWareing et~all ()2007l 
hereafter WZO) performed three-dimensional hydrody- 
namic simulations, from the beginning of the slow AGB 
wind phase to the end of the post-AGB/PN phase. They 
developed a "triple-wind" model including a slow AGB 
wind, a fast post-AGB wind, and a third wind reflecting 
the movement through the ISM. WZO found that the 
ISM interaction strongly affects the outer structures of 
PNs. The PN-ISM interaction process can be divided 
into four stages (WZO 1-4). In the first stage (WZO 

1) , the PN expands freely within the bubble of the undis- 
turbed slow AGB wind material and is unaffected by the 
AGB wind-ISM interaction. In the second stage (WZO 

2) , the PNs interact with the bow shock formed during 
the AGB wind-ISM interaction process, and the density 



and the emission in the upstream side increase accord- 
ingly. The third stage (WZO 3) is defined by the geo- 
metric center of the PN moving downstream away from 
the CS. In the fourth and final stage (WZO 4), the PN 
is completely disrupted and the CS is outside the PN. 

The relative velocity of the PN and the ISM density are 
important parameters for us to consider in order to model 
and understand the PN-ISM interaction processes. Both 
of these parameters are relevant to the Galactic orbits 
of PNs. The orbital parameters can also be used to dis- 
tingu ish different populations in the Galaxy t Wu et al.l 
I2009f ). but the available orbital parameters for PNs in 
the Galaxy are very rare. iRauch et al.l ()2004D calculated 
the orbit of PG 1034 -hOOl and found that t his PN is in an 
advanced stage of interaction with the ISM. lKerber et ahl 
([2004) calculated the orbits of four PNs and confirmed 
that the motion of the PN shell through the ISM is the 
origin of the PN-ISM interaction process. In this re- 
search, using improved observation data, we calculated 
the Galactic orbits for eight nearby PNs, which doubles 
the number of PNs whose orbits are available. Using the 
derived orbital parameters of those PNs, we also discuss 
the PN-ISM interaction processes in our sample in detail. 

We present the derived Galactic orbits of our PN sam- 
ple in Section 2. In Section 3, for each PN in our sample, 
using the derived kinematical data, we discuss the inter- 
action between the PN and the ISM. Our conclusions are 
given in Section 4. 

2. GALACTIC ORBITS 

2.1. The sample 

The position, distance, absolute proper motion, and 
radial velocity are essential data for deriving the Galac- 
tic orbit of a PN. First, we searched the literature to find 
PNs whose distance data are available. There are many 
indirect methods of PN distance determination. On the 
one hand, statistical distances can be underestimates. 
On the other hand, spectroscopic distances can be over- 
estimates (Benedict ct al. 200_9.). We consider only PNs 
whose distances a re deri v ed fro m direct parallax mea- 
surements. iHarris et al.l (j2007D presented trigonomet- 
ric parallaxes of the CSs of 16 nearby PNs. The me- 
dian error in their derived parallaxes is 0.42 mas which 
is smaller than th e ±1 mas t ypical error for most Hip- 
parcos parallaxes (|ES Al ITqQTI ) . Using the Hubble Space 
Telescope, the absolu te parallaxes of th e CSs of four 
PNs in t he sample of [ Harris d~an (|2007h were also de- 
rived bv iBenedict et al.l (I2009D. The derived errors m 

3 times smaller than those 



IBenedict et al 
in IHarris et al 



( 2009 ) are 2 
((2Ci07f) . 

Based on a lar ge number of niulti- wavelength data, 
iFrew et al.l (|2010t ) found that the nebula around the CS 
PHL 932 is a small H II region a nd not a real PN. Re ject- 
ing PHL 932 from the sample of IHarris et al.l (|2007[ ). we 
searched the proper motion and radial velocity data for 
the other 15 PNs in t he literature. Using a vailable all-sky 
astrometric catalogs, iKerber et al.l (|2008( ) collected reli- 
able proper motion information (with significance bet- 
ter than 3cr in at least one compon ent) for a sample o f 
234 PNs. Five PNs in the sample oflHarris et al.l (I2007D 
are not listed in the catalog of IKerber et al.l (|2008D : 
NGC 6853, RE 1738-H665, PG 1034-hOOl, A7, and A24. 
We used the VizieR catalog access tool to search the 
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proper motion data for these five PNs. No proper mo- 
tion data are available for NGC 6853, RE 1738+665, 
and A24. We fo und a proper motion for A7 in the 
UCAC 2 catalog ()Zacharias et al.ll2004[ ). but the errors 
are three times larger than the proper motions. We found 
a pro per motion for PG 10 34+001 in the UCAC 3 cat- 
alog ([Zacharias et al.|[2009D , but the proper motions are 
three times larger than their errors in both components. 
So, PG 1034+001 is included in the followi ng analysis. 

For the 10 PNs included in the catalog of lKerber et al.l 
(jlOOl), proper motions derived from different astromet- 
ric ca talogs are listed. Tho se catalogs inclu de Hippar- 
cos (ESAl ll997D Tvc ho 2 (IHog et al.l [2000l) UC AC 2 
(IZacharias^ eTaLl 120041) GSC 2 (iLas ker et alll2008l ). and 
USNO-B (|Monet et al.il2003D . The former four catalogs 
are all based on the International Celestial Reference Sys- 
tem (ICRS), but the proper motions of USNO-B were 
derived relative to the YS4.0 catalog and are not abso- 
lute. We do not adopt the data taken from the USNO-B 
catalog. 

The UCAC 3 catalog was made available in 2009. In 
addition to the early epoch data used for the UCAC 2 
proper motions, UCAC 3 utilizes data from about 5000 
unpublished astrographic plates as well as a complete re- 
reduction of the SPM first epoch data to significantly 
lower the systematic e rrors of UCAC 3 proper motions 
(jZacharias et al.ll2009[ ). The proper motions in UCAC 3 
are improved with respect to those in Hipparcos, Tycho 
2, and U CAC 2. Thus for CS with V < 14.0 mag in the 
sample of iHarris et al.l |2007,), the proper motions from 
UCAC 3 are adopted. For stars with V > 14.0 mag, 
there are larger system errors in the proper motions of 
U CAC 3. So, for CS with V > 14.0 mag in the sample 
of IHarris et al.l (|2007l ). the proper motions derived from 
the GSC 2 catalog and listed in Kerbcr et al. (2008) are 
adopted, but, the proper motions of NGC 6270 derived 
from UCAC 2 are very different from those derived from 
GSC 2, and no data can be found in the UCAC 3 catalog 
for this object. We rejected this PN in our analysis. 

The SIMBAD database was used to search the radial 
velocity data for the 10 PNs with known distances and 
proper motions. In the end, eight PNs with reliable 
distances, proper motions, and radial velocities were in- 
cluded in our final sample. In TableUthe posit ions (a, S) 
of each PN adopted from iKerber et al.l ( |2003| ) are listed 
in co lumns 2 and 3 , the V magnitudes of the CSs adopted 
from IHarris et all (f2007.) are listed i n column 4, c olumn 
5 li sts the derived parall axes from IHarris et al.l ()2007l ) 
and lBenedict et all (|2009[ ). proper motions (na cos 5, us) 
taken from the UCAC 3 and GSC 2 catalogs are listed 
in columns 6 and 7, the position angles (P.A.s) of proper 
motions are listed in column 8, and radial velocities and 
their references are listed in columns 9 and 10. 

In Table [2] the angular diameters of the eight 
PNs, except for PG 1034+001, are taken from 
iTweedv fc Kwitteii (|1996[) and are listed in column 2. 
Using the distances listed in column 4, the linear diam- 
eters of these PNs were derived and are listed in col- 
umn 3. The electron densities, JT-pji^, of PNs, if avail- 
able, are listed in column 5. The volume densities, 
""ISM' diffuse warm neutral medium (WNM) and 

warm ionized medium (WIM) at the observed positions 
of each PN were calculated according to the model of 



IKalberla fc Kerol (fT998h and are listed in columns 6 and 
7. The expansion velocities of each PN, if available are 
listed in column 8, and the corresponding references are 
listed in the final column. Using the linear diameter and 
expansion velocity, the kinematic age for each PN was 
calcul ated and is listed in column 9. IKalberla fc KerpI 
(119981 ) used a hydrostatic Galactic model, which assume 
that the interstellar magnetic field, gas, and cosmic rays 
form a dynamic equilibrium with the Galactic gravita- 
tional potential. They considered various Galactic gas 
phases in their model; these gases include cold neu- 
tral medium (CNM), WNM, WIM, neutral halo medium 
(NHM), and hot halo medium (HHM). Because the prob- 
ability of a PN-ISM interaction occurring for CNM is 
very small ()Borkowski et al.lll990D , and all of the PNs in 
our sample belong to the Galactic disk population (see 
the next section) , where the densities of NHM and HHM 
are much lower, these two gases need not be considered 
in the PN-ISM interaction for disk PNs. Hence we cal- 
culated the density for W NM and WIM only with the 
model of IKalberla fc Kerd ([1998). Table [2] indicates that 
the total densities, fijs]V[, of WNM and WIM around the 
eight PNs are close to ~ 0.1 cm~'^. 

2.2. Orbital parameters 

Using the data listed in Table [T] and adopting a so- 
lar motion (U. V W)r:, = (10.0, 5.2, 7.2) km s"! 
(jPehnen fc BinnevI 1998), the local standard of rest 
(LSR) velocities of PNs were determined. The LSR ve- 
locities were then corrected to the Galactic standard of 
rest (GSR) by adopting the galactocentric distance of the 
Sun Rq — 8.0 kpc and a rotation veloci ty of the LSR at 
the position of the Sun of 220 km s"! f Wu et al.l[200l 
and reference therein). The observed GSR positions and 
velocities of the eight nearby PNs are listed in Table [3l 
The position coordinates x, y, z and velocity components 
U, V, W refer to a galactocentric righthanded Carte- 
sian coordinate system with the directions pointing to 
the Galactic center, direction of rotation, and the Galac- 
tic north pole. The errors in the velocity components 
listed in Table [3] include the errors in distances, proper 
motions, and radial velocities of PNs listed in Table [TJ 
The LSR velocities of the PNs are listed in the last col- 
umn. 

We employed the axisymmetric Galactic gravitational 
potential models of | Allen fc Santillan (1991, hereafter 
AS91) and lFlvnn et al.l (|l996l hereafter FSC96). These 
two models consist of three components: bulge, disk, and 
spherical halo. Both models admit two conserved quan- 
tities: the total energy E and the z-component Jz of 
the angular momentum vector. The disks are v ery dif- 
ferent between the two models (jWu et al.ll2009| ). Using 
the initial conditions listed in Table |3l the orbits of the 
PNs were calculated backward in a time interval of 2 Gyr 
with the two given gravitat ional potent i al mo dels. The 
Bulirsch-Stoer algorithm of iPress et al.l ()1992f ) was used 
for the integration. The relative change in the total en- 
ergy over the 2 Gyr integration time is of the order of 
10-i4_ 

The orbital parameters e, zmax, and calculated 
with the AS91 and FSC96 models are hsted in Table 
m e = (i?a — Rp)/{Ra, + ^p) is the orbital eccentric- 
ity, where i?a and Rp are apogalactic and perigalactic 
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distances from the Galactic center, respectively, which 
are determined from the average maximum and min- 
imum galactocentric distances in the calculated orbits 
within the integration time of 2 Gyr. zniax is the aver- 
age maximum ve rtical distance above the Galactic plane 
(jWu et al.l[2?)n9h . is the mean time interval for the 
PN to cross the Galactic plan e from one Zma x to the 
other in the opposite direction (iWu et aLll2009D . The z- 
component of the angular momentum vector of each 
PNs is also listed in Table H For all eight PNs m our 
sample. Tables [2] and |4] indicate that the mean time 
for each PN is much larger than its kinematic age. 

The errors of orbital parameters listed in Table 2] were 
calculated based on a Monte Carlo method. Because the 
intrinsic uncertainties of the orbits within the 2 Gyr in- 
tegration interval are very small, only the errors in the 
input data were used to calculate the errors of the de- 
rived orbital parameters (jWu et al.l[2'009i ). For each PN, 
the initial condition was generated by adding Gaussian 
deviates to the observed proper motion, radial velocity, 
and distance. The errors for the input data listed in Ta- 
ble [3] are adopted as the standard deviations. For each 
PN, the errors for its orbital parameters were calculated 
based on 1000 separate integrations. 

The average difference for the orbital parameter e, 
-^max, and derived from the AS91 and FSC96 models 
are 0.006 ± 0.003, -54.1 ± 13.2 pc, and -10.5 ± 0.3 Myr, 
respectively. Table |4] shows that, the relative differences 
in the derived orbital parameters due to the adopted dif- 
ferent potential models are smaller than those derived 
from observational errors. The relative large differences 
in ^max and for the two models are due to the fact 
that the mass of the disk in t he FSC96 model is smaller 
than that in the AS91 model (|Wu et al.ll2"009| ). 

The orbits calculated with the AS91 and FSC96 models 
in the time interval of 2 Gyr are presented in Figure 
[U where p — \J +2/^- For all PNs in our sample, 
their orbits projected on the Galactic plane in Figure [T] 
clearly indicate the periodic motions of PNs in the plane. 
For most of the PNs, the meridional orbits in Figure [T] 
are of a boxy-like type and PNs move in the meridional 
plane within the limited regions almost filling the boxes. 
However, the meridional orbits of Sh 2-216 and PuWe 
1 in the AS91 model and that of PG1034-K001 in the 
FSC96 model indicate that they only move with special 
paths within the meridional plane. Thus, the weird orbits 
of Sh 2-216 and PuWe 1 derived from the AS91 model 
and that of PG1034-H001 derived from the FSC96 model 
are not used in the following analysis. These weird orbits 
are caused by the Galactic models adopted. Thus, we 
need to check each model for every object in detail, which 
is beyond the scope of this paper. This problem is under 
i nvestigation and will be discussed in a subsequent paper. 

iSchonrich et al.l (|2010( ) re-examined the velocity of the 
Sun with respect to the LSR. They found that the classi- 
cal determination of its component Vq in the direction of 
Galactic rotation is undermined by the metallicity gra- 
dient in the disk. Using a chemodynamical model which 
accounts for these effects, they obtained ([/, V , W)q — 
(11.1, 12.24, 7.25) km s~^ In particular, Vp. is 7 km s~i 
larger than that of iDehnen fc BinnevI ()1998l) adopted in 
this study. Using th e new parameters of the Sun derived 
bv iSchonrich et all (poTol ). we re-calculated the orbits of 



the eight PNs in our sample with the AS91 and FSC96 
models. Using the AS91 model, the mean differences for 
the orbital parameters e, Zrnax, a nd derived with the 
different parameters of the Sun (jSchonrich et al.l [20101 : 
IDehnen fc Binned [T998h are -0.018 ± 0.006, 3.9 ± 0.7 
pc, and 0.9 ± 0.1 Myr, respectively. The corresponding 
mean differences for e, zmax, and with the FSC96 
model are -0.018 ± 0.005, 2.0 ± 0.9 pc, and 1.5 ± 0.2 
Myr, respectively. The effects of adopting different pa- 
rameters of the Sun are very small for the derived orbital 
parameters with different Galactic models. 

2.3. Population membership 

iWu et afl (|2009[ ) analyzed the orbital eccentricity dis- 
tributions for different Galactic populations and found 
that no thin-disk star has an eccentricity greater than 
0.3 and no thick-disk star has an eccentricity less than 
0.1. The orbital eccentricities of the thin-disk and the 
thick-disk stars overlap between 0.1 and 0.3. The mean 
e for the thin-disk stars is ~ 0.1 and the mean for thick- 
disk stars is ~ 0.4. In our sample, only Sh 2-216 has 
e less than 0.1 which indicates it is a thin-disk PN. It 
is difficult to assign the population membership for the 
o ther seven PNs based on ly on their orbital parameter e. 

iPauli et all (|2003l 1200611 presented kinematic criteria 
for distinguishing thin-disk, thick-disk, and halo popula- 
tions in their DA white dwarf sample. Their kinematic 
criteria were deduced from a sample of F and G stars cho- 
sen from the literature, for which chemical criteria can 
be used to distinguish different populations. Their kine- 
matic population classification scheme is based on the 
position in the U — V velocity diagram, the position in 
the Jz ~ e diagram, and the Galactic orbits themselves. 
Among these three kinematic criteria, the Jz — e dia- 
gram can distinguish different populations more clearly 
(comparing Figures 2 and 3 of f Pauli et al. 2006)). In 
the Jz — e diagram of IPauli et al.i (|2006l . the thin-disk 
stars lie in a V-shaped area with e < 0.27 and Jz around 
1800 kpc km s~^. The thick-disk stars possess e > 0.27 
and Jz < 1500 kpc km s~^. Because of our adopted 
right-handed orientation of the coordinate system, the 
direct i on of Jz vector is opposite to that of IPauli et al.l 
(20031 12006D . We present the \Jz\ - e diagram in Fig- 
ure [21 the adopted orbital parameters e are derived from 
the AS91 model. Because of the weird orbits of Sh 2- 
216 and PuWe 1 in the AS91 model, their parameters 
e derived from the FSC96 model are adopted in Figure 
[21 The error of e for each PN in Figure [2l includes the 
error due to the uncertainty in the input data and that 
derived with different Galactic models. Due to the weird 
orbits of Sh 2-216, PuWe 1, and PG1034-I-001 in differ- 
ent models, their errors of e in different models cannot 
be confirmed. Thus, only errors derived from the input 
data are considered for the three above-mentioned PNs. 

In Figure [H dashed lines divide this figure into two 
regions for different populations. The top left region en- 
closed by these two dashed lines is the one for the thin- 
disk population. The other region is occupied by the 
thick-disk population. Figure [2 shows that PG 1034+001 
and A31 lie in the region of the thick-disk population and 
Sh 2-216, DeHt 5, NGC 7293, A21, and Ton 320 belong 
to the thin-disk population. PuWe 1 has low \Jz\ and 
e ~ 0.2 which indicates it likely belongs to the thick-disk 
population, but the large errors in its orbital parameters 
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make its population classification more uncertain. 
3. COMMENTS ON INDIVIDUAL OBJECTS 
3.1. Sh 2-216 

The orbit of S h 2-216 was first calculated by 
IKerber et all (I2004D with th e AS91 model . Usin g the 
method of lPauli et all (|2003[ ). IKerber et ahl (|2004D iden- 
tified Sh 2-216 as a thin-disk PN. They used the t rigono- 
metrical distance derived by [Harris et al.l (I1997D which 
is consistent with that derived by [Harris "etaiT (I2OOI . 
The proper motions (/ig , u s cos a ) = (21.7, —12.7) mas 
yr~^ used by IKerber et al.l ()2004[ ) were taken from the 
UCAC 2 catalog which are also consis tent with those 
taken from the UCAC 3 catalog. Ker ber et all (|2004l ) 
used the radial velocit y of 11.9 km s ~^ derived by 
iTweedv fc Napiwotzkil ll992D which is h alf of that de- 
rived bv iRauch et al.l (|2007D . iRauch et a l. (2007 ) have 
pointed out that the result of iTweedv fc Napiwotzkil 
(fl992 !) may be less accurate because the CS of this PN 
was possibly not well centered i n the observation aper- 
ture. Using the radial velocity of ITweedv fc Napiwotzkil 
(|1992D . we re-calculated the orbit of Sh 2-216 in the AS91 
model, the orbital parameters e, Zmax, and are de- 
rived as 0.03, 142 pc, and 33.3 Myr, respectively. The 
adopted different radial velocities of Sh 2-216 do not af- 
fect the derived orbital parameters with the same Galac- 
tic model. Using t he orbital parameters d e rived with the 
radial velocity of ITweedv fc Napiwotzkil (|1992D . based 
on its position in Figure O Sh 2-216 is still identified 
as a th in-disk PN, which is also consistent with that of 
IKerber et al. (200 4|). Therefo re, using the radial veloc- 
itv of pTweedv fc Napiwotzkil (jl992f) will not affect the 
population membership determination of Sh 2-216. 

The inte raction between Sh 2-216 and the ISM was 
studied by iTweedy et a.l (|1995[ ) in detail. A region of 
enhanced emission in the east of the PN and the dis- 
placement of the CS with respect to the center of the 
PN clearly indicate the PN-ISM interaction. The P.A. 
of the proper motions of Sh 2-216 also indicates that the 
direction of the movement of the PN is consistent with 
the position of the PN-ISM interaction region. Using the 
kinematic age 0.45 Myr and the W velocity of Sh 2-216 
listed in Table |3l this PN only moved 6.2 pc away from 
its birthplace in the vertical direction. We have pointed 
out in the previous section that the kinematic age of this 
PN is much smaller than the mean time for it to cross 
the Galactic plane. 

Assuming that the ISM in the vicinity of the Sun has 
the same rotation velocity of the LSR, the LSR velocities 
of the PNs listed in Table |3| can be considered as the ap- 
proximately relative velocities of those PNs with respect 
to their nearby ISM. Using a relative velocity of ~ 20 
km s~^, with an estimated electron density of njy^ ^ 
8 cm' ansom "eralll2008l) a nd an expansion velocity 
of ~ 4 km s-i ()Revnoldslll985l) . the models of BSS and 
SBS should need a local ISM density of njgjyj ~ 1.4 cm^'^ 
to explain the observed PN-ISM interaction of Sh 2-216. 
According to the model of Kalbcrla & Kcrp (1998), the 
density of the diffuse ISM at the position of Sh 2-216 is 
0.12 cm~'^ is much lower than that es t imate d from the 
models of BSS and SBS. ITweedv et a] ()1995[ ) estimated 
a relative velocity of only ~ 5 km s~^; an unusually high 
ISM density of r^jgjyj ^ 4 cm"'^ can be estimated with 



this much lower relative velocity. If we use the radial 
velocity of ITweedv fc Napiwotzkil (|1992| ). based on the 
models of BSS and SBS, the local ISM density can be 
derived as 2.2 cm""^, whic h is still much larger than that 
derived from the model of iKalberla fc KerpI ([1998'). 

Based on the above analysis, we find that the local 
density of ISM near Sh 2-216 derived from the mod- 
els of BSS and S BS is much larger than that derived 
fr om the model of IKalberla fc Kerpj (|1998f ). The model 
of lKalberla fc KerpT i 19981 ) represents the large-scale ISM 
distribution in the Galaxy, the values listed in Table [2] 
are the average and smoothed density of ISM at the po- 
sition of Sh 2-216. In fact, large d ensity perturbations 
on sm all scales of ISM are observed (iDickev fc LockmanI 
I1990D . hence the model of lKalberla fc Kerpj ()1998D mav 
underestimate the local density of ISM interacting with 
Sh 2-216. On the other hand, if the densi ty of ISM de- 
rived from the model of lKalberla fc K erpI (1 1998.) is cor- 
rect at the position of Sh 2-216, other factors must be 
considered in the models of BSS and SBS to explain the 
observ ed PN-ISM interaction for Sh 2-216. ITweedv et"a] 
(|1995l ) pointed out that the ISM magnetic field is needed 
to mold the interaction between Sh 2-216 and its nearby 
ISIvl 

iRansom et al.l (|2008[ ) presented 1420 MHz polarization 
images of Sh 2-216. An arc of low polarized intensity in 
the northeast portion was observed and was coincident 
with the optically identified PN-ISM interaction region. 
Seve ral polarization-angle "knots" appeared along the 
arc. IRansom et al.l ()2008| ) estimated the rotation mea- 
sure through the knots and derived a line-of-sight mag- 
netic field in the interaction region as 5.0 ± 2;0j£G;_This 
result is consistent with that estimated by ITweedv et a.l 
(IT995h . 

In the simulations of iVillaver et al.l (|2003l ) and WZO, 
with a low relative velocity ~ 20 km s"'^, the PN-ISM 
interaction can be observed for a lower ISM density of 
"■ISM ^ ^^-^^ cm~'^. The magnetic field is also not con- 
sidered in their simulations. The CS has been displaced 
from the center of the PN, so Sh 2-216 is in the WZO 3 
stage. 

3.2. PuWe 1 

Figure 2 of ITweedv fc Kwitted (|1996[ ) shows that 
PuWe 1 has a spherical enhanced emission shell interact- 
ing with the ISM, but the northwest sector of the shell is 
dimmer than the other portions. The P.A. of the proper 
motions of PuWe 1 indicates that this PN is moving to 
the southeast which is consistent with the direction of 
the PN-ISM interaction region. The orbit of PuWe 1 
shows that it is moving toward its maximum distance 
from the plane of the Galaxy. Using the kinematic age 
0.04 Myr and the W velocity of PuWe 1 listed in Table 
[3l this PN only moved 1.1 pc away from its birthplace in 
the vertical direction. 

The relative velocity of PuWe 1 is 50 km s^^ which 
is the typical velocity in the models of BSS and SBS 
for disk PNs. Using the ave rage volume fil ling factor of 
Galactic disk PNs e = 0.37 (jBoffi fc Stranghellinil ll994'). 
the mean elect ron dens i ty of PuWe 1 can be estimated 
as 3.1 cm ~ 3 (IPhillipslll998h . According to BSS and 
SBS, with an expansion velocity of 27 km s~^, the ISM 
density can be calculated as ~ 0.05 cm ~^ which is close 
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to the average diffuse ISM density njcj]^^ = 0.1 cm ^ 
derived from the model of iKalberla &: KerpI (|1998[ ) at 
the position of PuWe 1. Because the CS is still in the 
center of the PN, PuWe 1 is in the WZO 2 stage. 

3.3. A21 

Both the I P III] (jKwitter et al.l I1983D and Ha 

(|Hua fc Kwoki 11999^ images of A21 show that this PN 
is a one-side nebula, whose bright filaments are in the 
southeast and a faint rim is on the opposite side. The 
bright filamentary structure in the southeast indicates 
the PN-ISM interaction. The P.A. of the proper motions 
of this PN is consistent with the direction of the PN-ISM 
interaction regions. Tweedy & Kwittcr (1996) explained 
the one-side structure of A21 as the result of a highly 
inhomogeneous ISM in the vicinity of this PN. 

The orbits of A21 indicate that it is a typical thin-disk 
object. Using the derived kinematic age 0.02 Myr and a 
W velocity of A21, this PN only moved 0.2 pc away from 
its birthplace in the vertical direction. Different orbital 
parameters derived with different Galactic models cannot 
affect the PN-ISM interaction of A21. 

The electron densities of A21 in six regions of t he bright 
filaments were derived by IKwitter et al.l (|1983[ ) and the 
mean is 206 ± 42 cm~^. Using a relative velocity of ^ 30 
km and an expansion velocity of ~ 30 km s~^, ac- 
cording to BSS and SBS, the density of ISM in the vicin- 
ity of A21 c an be estimated a s ^ 5 cm~'^. According to 
the model of IKalberla fc Kerd (|1998[ ). the density of ISM 
at the position is 0.09 cm which is much smaller than 
that derived from the models of BSS and SBS. BSS has 
pointed out that an asymmetric H II region surround- 
ing A21 should exist, which would explain the observed 
shape of the P N and the estimate d high ISM density. The 
model of Kalberla fc KerpI (fT998h may underestimate the 
local density of ISM in the vicinity of A21. Because the 
CS of A21 has been displaced from the center, this PN 
is in the WZO 3 stage. 

3.4. Ton 320 

This PN was found during a sea rch around the hot 
H-rich white dwarf Ton 320 (|Tweedv & Kwittcr 199^. 
With a distance of 532 pc derived bv [Harris et al., (j2007| . 
this PN has a linear diameter of 4.6 pc, which is big- 

ger than Sh 2-216 . The [N II] image of Ton 320 

(jTweedv fc Kwitted [l99l shows that this PN is very 
faint and it has a similar one-side filamentary structure 
as that of A21. The bright filaments in the southeast 
indicate it is interacting with the ISM, no faint filaments 
can be observed in the opposite direction in the same 
[N II] image. The P.A. of the proper motions of Ton 320 
is consistent with the direction of the PN-ISM interaction 
region. 

Using the kinematic age 0.11 Myr and the W velocity 
of Ton 320, this PN only moved 1.1 pc away from its 
birthplace in the vertical direction. Because there is no 
electron density data available for Ton 320, we cannot 
estimate the ISM density with the models of BSS and 
BSB. However, using the diffuse ISM density ~ 0.1 cm~'^ 
at the position of Ton 320 derived from the model of 
IKalberla fc Kerd (|1998D . with a relative velocity of 40 
km and a mean expansion velocity of PNs ^ 20 km 
(|Weinberger.,1989.) . according to the models of BSS 



and SBS, the mean electron density of Ton 320 can be 
estimated as 3.6 cm^-^. This electron density is close to 
those of Sh 2-216 and PuWe 1. All three of those PNs 
have a large size and extremely low surface brightness, 
so the estimated very low electron density for Ton 320 
is reasonable. The very large size and very low density 
of Ton 320 also indicate that it has evolved to its later 
stage. The CS is displacing from the center of the PN, 
which indicates that Ton 320 is in the WZO 3 stage. 

3.5. A31 

Th e Ha and [N II] images of A31 (jTweedv fc KwitteH 
1 19941 ) show that its southern rim is brighter and sharp- 
edged, whereas the northern side is faint and diffuse, 
which indicates its interaction with the ISM. The orbits 
of A31 indicate it is a thick-disk object and is moving 
toward the Galactic plane. Using the kinematic age 0.13 
Myr and the W velocity of A31, this PN only moved 3.9 
pc away from its birthplace in the vertical direction. The 
P.A. of the proper motions of this PN is also consistent 
with the direction of the PN-ISM interaction region. 

Using the average filling factor of disk PNs e — 0.37 
(jBoffi fc Stranghelhnilfiggl . the mean electron density 
of A31 can be estimated as 11.4 cm'^ (jPhillipd [19981 ). 
Using the relative velocity of ^ 70 km s~^ and expan- 
sion velocity of ~ 10 km s^^, according to the models 
of BSS and SBS, the density of ISM in the vicinity of 
A31 can be estimated as 0.16 cm^'^ which is close to 
the aver age density ^tysM ^ 0.1 cm~'^ derived from the 
model of IKalberla fc Kerd (jT99l for the diffuse ISM at 
the position of A31. 

3.6. PG 1034+001 

The orbit of PC 1034-hOOl was calculated by 
IRauch et all (j2004f) with the AS91 Galactic model. The 
radial velo city and proper ri iotions of PG 1034-1-001 
adopted by IRauch et al.l (|2004D are consistent with those 
adopted in this paper, but the distance of 155 pc adopted 
by Rauch et al. (2004) is derived based on spectroscopic 
observation (| Werner et al.l 119951 ) and is less than our 
adopted di s tance of 210 pc. Using the method of 
iPauli et al] (j2003D and based on its position in the U 
- V v elocity diagram, - e diagram, and the Galactic 
orbit. IRauch et all (j200l found that PG 1034+OOljsa 
thin-disk star. Using the improved criteria of Pauli et al.l 
(1^006) for distinguishing different populations, we find 
that this star is a thick-disk object. The different pop- 
ulation membership determination is mainly due to the 
adopted different distances for this star. 

The largest known PN surrounding the DO white 
dwarf P G 1034-hOOl with a diameter of ^ 2° was first 
found bv lHewett et al.l (j2003[ ) based on the investigation 
of spectra from the Sloan Digital Sky Survey (SDSS; 
I York et alj I2000D . A careful inspection of the South- 
ern H - Alpha Sky Survey Atlas (SHASSA: lGaustad et all 
120011 ) has shown that this PN is surrounded by an el- 
liptical emission shell with a diameter of ~ 6°, and a 
fragment of a wider shell wi th a diameter of ^ 10° is 
also visible in the northeast (jRauch et al.ll200"^ . Using 
the kinematic age 0.54 Myr and the W velocity of PG 
1034+001 hsted in Tablejl this PN moved 10.2 pc away 
from its birthplace in the vertical direction; this distance 
is only half of its linear size. 
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iRauch et al.l (|2004f) noted that the major axis of the 
halo of the outer shell with a diameter of ~ 6° has a P.A. 
of ~ 310° which is consistent with the P.A. of the proper 
motions of 290°. They pointed out that this structure 
can be understood as an advanced PN-ISM interaction, 
but the P.A. of the proper motions is opposite to that of 
th e strongest emission r egion in the i nner nebula found 
bv lHewett eFall ()2003[ ). IRauch et all (|2004[) also pointed 
out that the DO white dwarf PG 1034+001 has probably 
experienced a final helium-shell flash in its past, which 
produces the second PN, but according to the simula- 
tions of WZO, the outer shell with a diameter of ~ 6° is 
probably the bow shock structure formed by the interac- 
tion between the slow AGB wind and the ISM; the inner 
shell with a diameter of ^ 2° is the PN shell formed by 
the interaction between the fast post- AGB wind and the 
slow AGB wind material. The PN is still within the bub- 
ble of undisturbed AGB wind material and unaffected by 
the ISM interaction, which indicates it is in the WZO 1 
stage. Because this PN is still in the WZO 1 stage and 
no electron density of this PN are available, the models 
of BSS and BSB cannot be used to explain the PN-ISM 
interaction and to derive the local density of ISM for this 
PN. 

3.7. NGC 7293 

NGC 7293 (Helix nebula) has a bright complex helical 
structure which is composed of at least two bipolar nebu- 
lae, with their poles at different orientations along P.A.s 
125° and 50°, respectively (jMeaburn et al.l |2008y. The 
very deep Ha-j - |N II] image of NGC 7293 obtained by 
iMeaburn et al.l (|2005f ) shows that there is a bow-shaped 
filament at the northeast of the bright helical filaments of 
this PN, and a faint jet-l ike feature at the op posite side 
of the PN is also visible. iKerber et al.l ()2004l ) calculated 
the orbit of NGC 7293 with the AS91 model and found 
that this PN is a thin-disk object. In this paper, we also 
identified NGC 7293 to be a thin-disk PN, w hich is con- 
sistent with the result of iKerber et al.l (|2004l ) . Using the 
kinematic age 0.03 Myr and the W velocity of PG NGC 
7293 listed in Table |3l this PN moved 0.4 pc away from 
it s birthplace in the v ertical direction. 

iTweedv fc Kwitteii (11996,) noted that the main shell 
of NGC 7293 does not show any signs of the PN-ISM 
interaction. They also pointed out that the asymmet- 
ric halo of this PN is due to the symmetrically ejected 
material from the CS which runs into the higher den- 
sity northeast region of its southwest counterpart. With 
the average electron densi ty of this PN ^ 630 cm"'^ 
(|Bofla fc Stranghellinilll994D . a relative velocity of 30 
km s~^, and an expansion velocity of ^ 20 km s~^, ac- 
cording to the models of BSS and SBS, the density of 
ISM in the vicinity of NGC 7293 can be estimated as 
^ 22 cm~^ which is mu ch larger than nyc;]\/|- = 0.08 de- 
rived from the model of lKalberla fc: Kerpl (jl998f ) for the 
diffuse ISM. 

According to the simulations of WZO, the bow-shaped 
structure in the halo is probably formed by the interac- 
tion between the AGB wind and the ISM. The P.A. of 
the proper motions of the PN is consistent with the outer 
bow-shaped region, which indicates a probable wind-ISM 
interaction. The inner bright helical shell is not affected 
by the ISM interaction, and this PN is in the WZO 1 



stage. 

3.8. DeHt 5 

The Ha and [N II] images of DeHt 5 

(jTweedv &: Kwitteii I1996D show a large amount of 
diffuse material surrounding the bright filaments of this 
PN, which indicates that the PN-ISM interaction has 
destroyed the intrinsic morphology of this PN. The 
very low expansion velocity of 5 km s^"'^ also indicates 
that DeHt 5 has evolved into its later stage. Using the 
kinematic age of 0.09 Myr and the W velocity of DeHt 
5 listed in Table [Sj this PN moved 1.0 pc away from its 
birthplace in the vertical direction. 

Using the relative velocity of ^ 60 km s^^ and the 
expansion velocity of 5 km s~^, according to the models 
of BSS and SBS, and using the diffuse ISM density ^ 1 
cm~^ derived from the model of lKalberla &: KerpI ()1998l) , 
the mean electron density of DeHt 5 can be estimated as 
4.5 cm~^ which is close to those of Sh 2-216, PuWe 1, 
and Ton 320. This very low density also indicates that 
DeHt 5 has evolved into its later stage. 

4. DISCUSSION AND CONCLUSIONS 

The orbits of eight nearby PNs are calculated with the 
AS91 and FSC96 Galactic models. The small number 
of PNs in the present sample is limited by the avail- 
able observation data for PNs in the Galaxy. In our 
present sample, only PNs with direct parallax distances 
are considered. The errors of the derived orbital parame- 
ters are calculated with a Monte Carlo method. We also 
compared the orbital parameters derived with different 
Galactic models. We find that the uncertainty of the de- 
rived orbital parameters is mainly from the observational 
errors in the input data. 

According to the kinematic criteria of IPauli et al.l 
(|2003l I2006D for distinguishing different populations in 
the Galaxy and using the — e diagram, we find that 
Sh 2-216, DeHt 5, NGC 7293, A21, and Ton 320 belong 
to the thin-disk population and PG 1034+001 and A31 
belong to the thick-disk population. PuWe 1 probably 
belongs to the thick-disk population, but with the large 
errors of its derived orbital parameters, its population 
classification becomes very uncertain. The present sam- 
ple of PNs is very small, hence we cannot come to any 
conclusion about the population distribution of Galactic 
PNs in the disk. Population classification for the eight 
PNs in our sample indicates that most of the PNs formed 
are closer to the Galactic plane. The dramatic enlarge- 
ment of the sample of PNs with reliable parallaxes will 
only be expected when the data of the Gaia mission are 
available. The primary objective of the Gaia mission is 
to survey more than one billion stars between 6 and 20 
mag. The expected accuracies for parallaxes, proper mo- 
tions and positions are 7 - 2 5 iias dowii to 15 mag and a 
few hundred /xas at 20 mag (lLindegrenll2010D. The mag: - 
nitudes of 96% CS in the catalog of IKerber et al.l ()2003h 
are less than 20.0; the precise distances and proper mo- 
tions for those stars obtained with the Gaia mission can 
be used to study their kinematics and dynamics in the 
Galaxy. 

We analyzed the interaction between the PNs in our 
sample and their nearby ISM. The PN-SIM interaction 
are observed in the eight PNs in our present sample. The 
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P.A.s of the proper motions of the PNs are also consis- 
tent with the directions of the PN-ISM interaction re- 
gions. For ah eight of the PNs in our sample, their kine- 
matic ages are much smaller than the time Tz for them 
to cross the Galactic plane, and the distances that they 
moved from their birthplaces are all smaller than their 
linear size. Different orbits derived with different Galac- 
tic models cannot affect the PN-ISM interaction for each 
PN. 

There are two main models to explain the PN-ISM in- 
teraction. One is presented by BSS and SBS. In their 
models, the PN interacts with the ISM directly, and the 
PN-ISM interaction can be observed only when the den- 
sity of the PN drops below some critical value. The mag- 
netic field is also important in the models of BSS and 
SBS. The other one is the triple-wind" model presented 
by WZO. In their model, the slow AGB wind formed in 
the early evolution phase of the AGB star interacts with 
the ISM, and then the PN interacts with the bow shock 
formed during the AGB wind-ISM interaction. So, in the 
model of WZO, the PN-ISM interaction can be observed 
in the early evolution phase of the PN. The magnetic 
field is not considered in the model of WZO. For both 
models, the relative velocity of the PN and the density 
of the ISM in the vicinity of the PN are very important 
parameters for studying the PN-ISM interaction. 

We use the above-mentioned models to explain the 
observed PN-ISM interaction in the eight PNs in our 
sample. Using the models of BSS and SBS, the local 
densities of ISM in the vicinity of PuWe 1 and A31 are 
derived and are consistent with those derived from the 
model of lKalberla fc Kerd ()1998D . For Ton 320 and DeHt 
5, no electron density data are available. Using the lo- 
cal densities o f ISM for these two PNs derived from the 
model of Kalb erla fc Ke rp (1998), the electron densities 
of these two PNs are derived from the models of BSS 
and SBS. Compared with the observed electron densi- 
ties of other PNs which have evolved into their later 
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stages, the derived electron densities of these two PNs 
are reasonable. For the above four PNs, the models of 
BSS and SBS yie ld consistent results from the model of 
Kalb erla fc KerpI ([1998), but for Sh 2-216 and A21, the 
derived local densities of ISM in the vicinity of these two 
PNs with the models of BSS and SBS are larg er than 
those derived from the model of lKalberla fc Kerg ([1998). 
For Sh 2-216, considering a magnetic field in the position 
of this PN in the models of BSS and SBS can explain the 
observed PN-ISM interaction. For A21, an asymmetric 
H II region surrounding this PN can explain the pre- 
di cted high ISM density. F or these two PNs, the model 
of lKalberla fc KerpI (|1998D will underestimate the loca l 
densities of ISM. The model of lKalberla fc Kerd ([l998l ) 
represents the large-scale and average ISM distribution 
in the Galaxy; it cannot predict the local density of ISM 
very precisely. Therefore, with the models of BSS and 
SBS, the local ISM densities in the vicinity of evolved 
PNs can be derived with th e observed PN-ISM interac- 
tion ()Borkowski et al.lll990f ). Because PG 1034-^001 and 
NGC 7293 are still in their early stage of evolution, the 
models of BSS and SBS can not used to derive the ISM 
densities for these two PNs. On the other hand, accord- 
ing to the model of WZO, Sh 2-216, A21, and Ton 320 
are in the WZO 3 stage, PG 1034+001 and NGC 7293 
are in the WZO 1 stage, and PuWe 1 is in the WZO 2 
stage. 
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Fk;. 1. Meridional Galactic orbits and orbits projected on to the Galactic plane in the time interval of 2 Gyr for the eight PNs in our 
sample calculated with the AS91 and FSC96 models. 
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TABLE 1 

The observed data of eight nearby PNs in our sample. 



Name 


a (2000.0) 


<5 (2000.0) 


V 


'^abs 


Ha COS (5 




P.A. 


^rad 


Reference 




(hh 


: mm : ss) 


(dd ; 


mm 


; ss) 


(mag) 


(mas) 


(mas yr"-"-) 


(mas yr"-"-) 


(°) 


(km s ) 




Sh 2-216 


04 


43 


21.3 


+46 


42 


05.8 


12.630 


7.76 + 0.33 


22.9 + 1.8 


-11.6 + 3.5'" 


116 


20.6 + 1.5 


7 


PuWe 1 


06 


19 


34.3 


+55 


36 


42.3 


15.534 


2.74 + 0.31 


19.6 + 14.0 


-24.8 + 8.0 


141 


18.3+ 1.2 


3 


A21 


07 


29 


02.7 


+13 


14 


48.8 


15.962 


1.85 ± 0.51 


2.4 + 2.0 


-10.0 + 2.0 


166 


28.8 + 5.2 


2 


Ton 320 


08 


27 


05.6 


+31 


30 


08.5 


15.702 


1.88 + 0.33 


-4.5 + 2.0 


-13.6 + 2.0 


198 


33.8 + 3.5 


4 


A31 


08 


54 


13.2 


+08 


53 


53.0 


15.519 


1.61 + 0.21=^ 


-18.4 + 2.0 


-16.5 + 2.0 


228 


44.0 + 1.0 


5 


PG 1034+001 


10 


37 


03.8 


-00 


08 


19.2 


13.211 


4.75 + 0.53 


-82.4 + 10.1 


31.2 + 10. 0*= 


290 


50.8+ 1.3 


6 


NGC 7293 


22 


29 


38.5 


-20 


50 


13.7 


13.525 


4.64 + 0.27=" 


27.4 + 1.4 


-5.4+1.5'' 


101 


-28.2 + 3.0 


2 


DeHt 5 


22 


19 


33.7 


+70 


56 


03.3 


15.474 


2.90 + 0.15=^ 


-17.4 + 2.0 


-19.8 + 2.0 


221 


-35.4 + 2.5 


1 



Referen ces. — (1) 'Barstow et al.' f^OOll); (2) IDurand et al.l ^\m^): (3 ') IGieseking eTal] lITgM) : (4) IGood et al.1 ||2005| '): (5) 
Iffippclcin fc Weinberge r (1990); (6) Holbcrg et at] II1998I ): (7) IRauch et al.l II2007I) 

Taken from Table 5 of iBenedict et al.l 1 1201)91) . 
Taken from the UCAC 3 catalog. 



TABLE 2 

The sizes, electron densities, and expansion velocities of eight nearby PNs in our sample and the number densities of ism at the 

positions of these PNs. 



Name 


angular diameter 
(arcmin) 


linear diameter 
(pc) 


distance 
(pc) 


(cm ^) 


njsM (WNM) 
(cm ^) 


njgM (WIM) 

(cm '^) 


tJexp 
(km s^^) 


kinematic age 
Myr 


Reference 


Sh 2-216 


100 


3.7 


128.9 


8.0 


0.10 


0.02 


4 


0.45 


5 


PuWe 1 


20 


2.1 


365.0 


3.1 


0.08 


0.02 


27 


0.04 


1 


A21 


9 


1.4 


540.5 


206.0 


0.07 


0.02 


32 


0.02 


4 


Ton 320 


30 


4.6 


531.9 




0.05 


0.02 




O.ll*" 




A31 


16 


2.9 


621.1 


11.4 


0.04 


0.02 


11 


0.13 


2 


PG 1034+001 


360^ 


22.0 


210.5 




0.07 


0.02 




0.54'" 




NGC 7293 


18 


1.1 


215.5 


630.0 


0.06 


0.02 


21 


0.03 


3 


DeHt 5 


9 


0.9 


344.8 




0.08 


0.02 


5 


0.09 


1 



References. — (1) IGieseking et al.l itTosg i : f 21 IHippelein fc Weinberger l(T990l '): (3') IMeaburn et al.ll|2008h : (41 IRecillas-Cruz fc Pisiiria l[T98lh : (5) 
IRevnoldsl ijigsl ) 

Taken from IRauch et al.l l[200l ). 

Use the mean expansion velocity of Galactic PNs 20 km s"! IIWeinbergeil[T989t) . 



TABLE 3 

The present positions and velocities of eight nearby PNs in the Galaxy. 



Name 


X 


y 


z 


U 


y 


W 


"LSR 




(pc) 


(pc) 


(pc) 


(km s"-') 


(km s"-') 


(km s"'^) 


(km s ) 


Sh 2-216 


-8119.9 + 5.1 


47.3 + 2.0 


1.1 + 0.0 


-14.4 + 1.6 


219.2 + 1.8 


13.3 + 1.6 


19.6 + 2.4 


PuWc 1 


-8324.2 + 36.7 


124.9+ 14.1 


111.9 + 12.7 


-20.8 + 7.1 


180.8 + 16.8 


27.3 + 22.0 


52.1 + 13.4 


A21 


-8474.3 + 130.7 


-222.6 + 61.4 


133.0 + 36.7 


-5.5 + 5.8 


189.4 + 8.4 


9.2 + 5.3 


32.4 + 3.3 


Ton 320 


-8436.9 + 76.7 


-88.1 + 15.5 


290.3 + 51.0 


-22.0 + 4.1 


187.5 + 7.5 


9.6 + 5.4 


40.4 + 3.5 


A31 


-8411.7 + 53.7 


-335.0 + 43.7 


322.6 + 42.1 


-35.2 + 4.9 


163.4 + 7.0 


-29.9 + 9.3 


73.1 + 5.7 


PG 1034+001 


-8054.0 + 6.0 


-130.8 + 14.6 


155.8 + 17.4 


-86.9 + 13.5 


197.0 + 7.9 


18.6 + 7.4 


91.8 + 12.9 


NGC 7293 


-7905.5 + 5.5 


69.0 + 4.0 


-181.0 + 10.5 


-23.6 + 2.2 


204.9 + 1.9 


15.5 + 2.8 


32.0 + 2.1 


DeHt 5 


-8121.5 + 6.3 


315.1 + 16.3 


69.6 + 3.6 


60.6 + 3.7 


209.9 + 2.8 


-10.7 + 3.3 


62.4 + 3.7 
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TABLE 4 

The Galactic orbital parameters of eight nearby PNs in our sample calculated with the AS91 and 

FSC96 MODELS. 



Name 






e 




(pc) 


T^Myr) 




(kpc km s~^) 


(AS91) 


(FSC96) 


(AS91) 


(FSC96) 


(AS91) 


(FSC96) 


Sh 2-216 


-1779.2 ± 14.5 


0.05 ±0.01 


0.05 ±0.01 


143.0 ±20.3 


190.0 ± 27.1 


33.9 


44.4 


PuWe 1 


-1502.4 ± 139.0 


0.19 ±0.07 


0.18 ±0.06 


335.0 ±321.3 


439.0 ±441.1 


35.1 


46.5 


A21 


-1606.3 ± 55.2 


0.14 ±0.04 


0.14 ±0.04 


158.0 ±43.7 


180.0 ± 58.2 


31.7 


41.7 


Ton 320 


-1583.9 ± 54.1 


0.16 ±0.03 


0.16 ±0.03 


287.0 ±38.5 


299.0 ±44.0 


35.1 


44.4 


A31 


-1386.3 ± 51.5 


0.26 ±0.03 


0.25 ±0.02 


504.0 ± 154.6 


620.0 ±212.6 


38.5 


50.0 


PG 1034+001 


-1598.0 ±64.2 


0.28 ± 0.04 


0.26 ± 0.03 


254.0 ± 72.5 


311.0 ± 101.5 


35.7 


46.5 


NGC 7293 


-1618.2 ±15.2 


0.10 ±0.01 


0.10 ±0.01 


248.0 ± 22.7 


290.0 ± 33.5 


34.5 


44.4 


DeHt 5 


-1723.8 ± 23.9 


0.17 ±0.01 


0.16 ±0.01 


130.0 ± 31.2 


163.0 ± 43.6 


33.9 


44.4 



1800 



1700 



1500 



4^ 



PG 1034+001 



0.00 



0.05 



0.10 



0.15 0.20 
e (AS91) 



0.25 



0.30 



0.35 



Fig. 2. — \ Jz\ — e diagram of the eight PNs in our sample. The top left region enclosed by two dashed lines is occupied by the thin-disk 
population. 
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